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ABSTRACT 

We report on the characterization and target ana- 
lysis of the small (s)RNA 162 in the methanoarchaeon 
Methanosarcina mazei. Using a combination of 
genetic approaches, transcriptome analysis and 
computational predictions, the bicistronic 
MM2441-MM2440 mRNA encoding the transcription 
factor MM2441 and a protein of unknown function 
was identified as a potential target of this sRNA, 
which due to processing accumulates as three 
stabile 5' fragments in late exponential growth. 
Mobility shift assays using various mutants verified 
that the non-structured single-stranded linker 
region of sRNA 162 (SLR) base-pairs with the 
MM2440-MM2441 mRNA internally, thereby 
masking the predicted ribosome binding site of 
MM2441. This most likely leads to translational re- 
pression of the second cistron resulting in 
dis-coordinated operon expression. Analysis of 
mutant RNAs in vivo confirmed that the SLR of 
sRNAi 62 is crucial for target interactions. 
Furthermore, our results indicate that sRNA 162 - 
controlled MM2441 is involved in regulating the 
metabolic switch between the carbon sources 
methanol and methylamine. Moreover, biochemical 
studies demonstrated that the 5' end of sRNA 162 
targets the 5-untranslated region of the c/s- 
encoded MM2442 mRNA. Overall, this first study of 
archaeal sRNA/mRNA-target interactions unraveled 
that sRNA 162 acts as an antisense (as)RNA on cis- 
and frans-encoded mRNAs via two distinct 
domains, indicating that c/s-encoded asRNAs can 
have larger target regulons than previously 
anticipated. 



INTRODUCTION 

In recent years, an increasing number of non-coding 
RNAs have been shown to participate in various regula- 
tory cellular processes in both pro- and eukaryotes, mainly 
acting as post-transcriptional ribo-regulators. Although 
the abundant class of eukaryotic miRNAs mainly act by 
base-pairing to the 3'- untranslated region (UTR) or 
coding sequence (CDS) of the cognate target-mRNA, 
most of the functionally characterized bacterial small 
regulatory RNAs (sRNAs) target the 5'-UTR of 
mRNAs (1-3). Besides several sRNAs which directly 
modulate protein activity, the majority of functionally 
characterized sRNAs from bacteria belong to the class 
of trans-encoded antisense RNAs. In most cases, they 
bind with only partial sequence complementarities to the 
5'-UTR of their target genes, and lead to translational 
repression by masking the ribosome binding site (RBS). 
Consequently, the association of the ribosomal 30S 
subunit to the RBS is inhibited, which is often coupled 
to target destabilization by RNases (2). However, repres- 
sive translational control and induction of mRNA degrad- 
ation is not limited to direct base-pairing with the RBS, 
because some sRNA have been shown to target sequences 
that are located far upstream of the RBS, within the CDS 
or in the intergenic region of polycistronic mRNAs (4-9). 
Besides repression of target genes, bacterial sRNAs have 
also been demonstrated to up-regulate gene expression by 
disruption of inhibitory secondary structures which se- 
quester the RBS, known as anti-antisense mechanism 
(10-12). The underlying sRNA-target interactions 
together with the coupled sRNA decay are often facili- 
tated by the bacterial RNA chaperone Hfq, which 
contains a Sm-like domain (13-15). Additionally, bacterial 
sRNAs have been shown to target multiple, often func- 
tionally related targets, some of which are regulators 
themselves, and a single sRNA can encompass the expres- 
sion of large regulons (reviewed in 16,17). 
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Although first examples of so-called m-encoded anti- 
sense (as)RNAs, which are encoded on the opposite strand 
of their target genes, were studied already in the early 
1980s on plasmids and transposons (18-20), asRNAs are 
less well characterized in comparison with their trans- 
acting counterparts. Such asRNAs generally act on their 
targets independently of Hfq, possibly due to the high 
extent of base-pairing between the sRNA and their 
target mRNAs. They overlap with their targets either at 
the 5'-UTR or 3'-UTR or within the center of the target 
mRNA. For asRNAs diverse regulatory mechanisms have 
been described, e.g. alteration of target stability, interplay 
with RNases, or interference with transcription (for details 
see 21 and 22). The massive detection of antisense tran- 
scription in various genome-wide transcriptome studies 
clearly suggests a more fundamental role in prokaryotic 
cells; e.g. in Helicobacter pylori for 46% of all ORFs anti- 
sense transcription was detected (23). Although a massive 
antisense transcription was also demonstrated in other 
human pathogens such as Listeria monocytogenis or 
Staphylococcus aureus, and also in the cyanobacterium 
Synechocystis sp. PCC 6803, the total impact of asRNA 
regulation is still not yet understood (24-27). 

In contrast to bacteria until recently no sRNAs had been 
identified in archaea, disregarding the role of the extensively 
studied eukaryotic like small nuclear RNAs (snoRNAs), 
which participate in ribosomal RNA biogenesis and 
tRNA maturation (28-31). Currently, several archaea 
have been examined for the presence of sRNAs on a 
genome-wide scale either by RNomics approaches or 
using high throughput sequencing (HTPS) techniques, e.g. 
RNA-Seq approaches (HTPS of cDNAs), resulting in the 
detection of high numbers of sRNA candidates (32-35). To 
study the impact of sRNAs and to get a deeper insight into 
transcriptional and post-transcriptional regulation in the 
methanogenic archaeon Methanosarcina mazei G61, specif- 
ically in response to nitrogen (36), we have recently applied 
a newly developed differential RNA-Seq (dRNA-Seq) 
approach for the analysis of primary transcriptomes (23). 
This approach resulted in the identification 248 sRNA can- 
didates, a high number of which were confirmed by 
northern blot analysis (36). Here, we report on the func- 
tional characterization and target analysis of the abundant 
sRNA 162 . Using genetic, biochemical and computational 
approaches, we demonstrate that sRNA 162 interacts with 
both, a cis- and a /ra«.s-encoded target mRNAs via two 
distinct domains, a mechanism which has not been shown 
for any studied sRNA in prokaryotes so far. 

MATERIALS AND METHODS 

Strains and plasmids 

Strains and plasmids, which were used in this study, are 
listed in Supplementary Table SI. Plasmid DNA was 
transformed M. mazei as described by Ehlers et at. (37). 

Construction of M. mazei mutants and generation 
of plasmids 

pRS699 was constructed by polymerase chain reaction 
(PCR) amplification of sRNA 162 including its native 



promoter and terminator from genomic M. mazei DNA 
using primers sl62-XhoI.for and sl62-KpnI.rev (Supple- 
mentary Table S2). The 628 bp PCR-fragment was TA 
cloned into the pCR4-TOPO (Invitrogen Karlsruhe, 
Germany), resulting in pRS699. After restriction with 
Xhol and Kpnl the fragment was inserted into the 
multiple cloning site of pWM321 (38). The resulting 
plasmid (pRS474) and mutant derivatives (discussed 
later) were transformed into M. mazei by liposome- 
mediated transformation as previously described (37,39). 
Puromycin-resistant transformants were selected as 
colonies that grew on minimal medium plates with 
trimethylamine as the carbon and energy source plus 
5 ug puromycin/ml during incubation. 

Using pRS699 as template the sRNA 162 Ml and M2 
mutants were generated by site-directed mutagenesis 
with the primer sets sl62-Mutl -for/rev or sl62-Mut2- 
for/rev resulting in pRS702 and pRS704. Further 
mutant derivatives were constructed by inverse PCR 
using pRS699 as PCR template: Replacement of the 
linker region between SL2 and SL3 by poly(T) with 
primers sl62-Mut3-for/rev (pRS706) and the deletion of 
the 3' end retaining the 5' fragment (nucleotide 1-65, 
pRS701) and 3' fragment (nucleotide 65-191, pRS708) 
of sRNA 162 with sl62-Mut4-for/rev and sl62-Mut5-for/ 
rev, respectively. The amplified products were religated 
and cloned into E. coli. Using this cloning strategy, the 
cognate sRNA variants kept their native sRNA 162 
promoter and terminator sequences. By using the afore- 
mentioned Xhol and Kpnl restriction sites, the mutated 
sRNA 162 variants were inserted into pWM321 as 
described earlier (pRS701, pRS703, pRS705, pRS707, 
pRS709). 

For the construction of the chromosomal sRNA 162 , 
deletion mutant by homologous recombination in 
M. mazei, plasmid pRS650 was constructed as follows. 
The flanking upstream region (~1 kbp) of sRNA 162 was 
amplified by PCR using sl62 Dell Xhol and sl62 Del2 
EcoRF The flanking downstream region (~1 kbp) was 
amplified by PCR using si 62 Del3 EcoRI and si 62 Del4 
Xbal. Both PCR products were restricted with the 
indicated restriction enzymes and separately ligated to 
correspondingly linearized pBluescript SK+ resulting in 
pRS648 and pRS649, respectively. The ~1.8 kbp EcoRI 
fragment of pRS207 carrying the pac cassette from 
Methanococcus voltae (40) conferring puromycin resist- 
ance was ligated in a three body ligation together with 
the upstream fragment (XhoI/EcoRI restricted) into the 
XhoI/EcoRI restricted plasmid pRS649 generating 
plasmid pRS650. pRS650 was transformed into 
M. mazei using a liposome-mediated transformation 
protocol as described previously (37). Southern blot 
analyses of genomic DNA from puromycin-resistant 
transformants were used to verify pac insertion as 
described by Ehlers et al. (37). 

The T7-hammerhead ribozyme fusion of sRNA 162 and 
sRNAi 6 2A63-88 was generated by commercial gene synthe- 
sis (Eurofins MWG Operon, Ebersberg, Germany). TA 
cloning of both of the products into vector pCR4-TOPO 
(Invitrogen Karlsruhe, Germany) yielded plasmids 
pRS765 and pRS766, respectively. 
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To fuse a T7 promoter to MM2440-41, the operon was 
amplified with genomic DNA as template as described (see 
in vitro transcription). The PCR product was TA cloned 
into pCR4-TOPO (pRS767) and the obtained plasmid was 
subjected to site-directed mutagenesis with primers 2441- 
com_Mutl for/rev (pRS768) and 2441-com_Mut2for/rev 
(pRS769) to obtain compensatory mutants. 

RNA isolation 

Exponentially or stationary phase cultures (50-70 ml) were 
harvested at 4°C and RNA was isolated by phenol extrac- 
tion as recently described (36) or using the RNeasy Midi 
Kit according to the manufacturer's instructions (Qiagen, 
Hilden, Germany) and published for M. mazei (41,42). 

Northern blot analysis 

Cells were grown under different nitrogen (N) 
availabilities and harvested at different growth stages (N 
sufficiency: early exponential phase, OD 60 o = 0.2; mid- 
exponential phase, OD 600 = 0.5, stationary phase, 
OD 60 o = 0.7; N limitation: mid-exponential phase, 
OD 60 o = 0.2). Northern blot analysis followed using the 
recently described protocol (36). sRNA 162 and its poten- 
tial processed derivatives were detected with a 5'-end 
radioactive labeled ssDNA oligo probe (Supplementary 
Table S2) and 5S rRNA as described in Jager et al. (36). 
Full length riboprobe preparation followed the same pro- 
cedure as described for T7 in vitro transcription, except 
that the MAXIscript T7 kit (Ambion; Applied 
Biosystems, Foster City, CA) was used according to the 
manufacturer. To amplify a T7 promoter containing PCR 
product which serve as template for antisense probe syn- 
thesis the primer pairs asPsl62 for and asPsl62rev, or 
asPsl71 for and asPsl71rev, were used, respectively 
(Supplementary Table S2). 

Rapid amplification of cDNA ends analysis 

RACE (rapid amplification of cDNA ends) analysis was 
performed to determine the transcriptional start site (TSS) 
as well as the transcript termination site of sRNA 162 . The 
5'-RACE system (Invitrogen, Karlsruhe, Germany) was 
used according to manufacturer's instructions with oligo- 
nucleotide 5'-RACE 2441 as gene-specific primer to deter- 
mine the TSS of sRNA 162 (Supplementary Table S2). 
3'-RACE analysis followed the 3'-RACE System 
(Invitrogen, Karlsruhe, Germany) according to manufac- 
turer's instructions after polyadenylation of RNA with 
Poly(A)-polymerase (NEB, Schwalbach, Germany). The 
precipitated RNA was resuspended in 1 1 -ul DEPC-H 2 0 
and directly used for 3' RACE using the oligonucleotide 
3'-RACE 2441 #2 as gene-specific primer (Supplementary 
Table S2). PCR products were subjected to TOPO-TA 
cloning using the TOPO-TA Cloning Kit (Invitrogen, 
Karlsruhe, Germany). 5' end and 3' end of sRNA 162 
were determined by DNA sequencing of both strands. 
For 5'-end analysis of the operon MM2440-41, the 
FirstChoice RLM-RACE kit (Ambion; Applied 
Biosystems) was used, according the manufacturer's 
instructions. The oligonucleotides 5'-2441-RLM (i?NA 
/inker «?ediated)-out and 5'-2441-RLM-in were used as 



gene-specific primers to determine the TSS. In the same 
way, the oligonucleotides 5'-2442-RLM-out and 5'-2442- 
RLM-in or 5'-2446-RLM-out and 5'-2446-RLM-in were 
used to determine the TSS of MM2442 or MM2446, re- 
spectively (Supplementary Table S2). 

In vitro T7 transcription, purification and 5' -end 
labeling of RNA 

Templates for in vitro transcription were either amplified 
from genomic M. mazei wild-type (wt) DNA or from 
plasmids carrying the respective constructs (Supple- 
mentary Table SI). Primer sequences are summarized in 
Supplementary Table S2. Plasmids carrying sRNA 162 
(or derivatives, discussed earlier) were PCR amplified sim- 
ultaneously fusing a T7 promoter to the PCR product. 
For the 3' end of sRNA 162 , the primer pair 
T7-sl62-SHORT-for and sRNA162-T7.rev was used, 
whereas the other T7 sRNA 162 fusions were amplified 
with T7-HH-sl62-for and sRNA162-T7.rev instead. 
Similar, T7-2440-41-for and T7-2440-41-revl or 
T7-2440-41-rev2 were used for the T7 fusions of the long 
and short versions of mRNA MM2440-2441, respectively. 
T3-MM2446-for and rev were used for the amplification 
and the incorporation of a T3 promoter to MM2446. 
In vitro transcription was performed using the 
MEGAscript T7 (or T3) kit (Ambion; Applied 
Biosystems) according to the manufacturer. Following ex- 
traction with phenol:chloroform:isopropanol (25:24:1 v/v) 
and purification with G-25 column, the RNA was 
precipitated overnight at — 20°C with 3 vol. of ethanol 
and 20 ug glycogen (Applied Biosystems). RNA quality 
and integrity were checked on a denaturing polyacryl- 
amide gel (6% PAA, 7M Urea). The in vitro transcribed 
RNA was dephosphorylated with FastAP (Thermosen- 
sitive Alkaline Phosphatase, MWG Fermentas), according 
to the manufacturer and radioactively labeled at the 5' end 
as described in Jager et al. (36), additionally supplement- 
ing the labeling reaction with 40 U of RNasin (Promega, 
Mannheim, Germany). 

A hammerhead ribozyme transcriptional fusion 
with sRNAi 62 

Usually, in vitro transcripts were created by incorporating 
a T7 promotor into the oligonucleotides used for PCR of 
the template or by cloning of the respective PCR product 
into a T7 promotor containing vector. To ensure efficient 
transcription with high yield, the PCR-based T7 fusion 
usually introduces 3 Gs at the TSS (+1 site). To prevent 
the addition of non-sequence specific nucleotides into the 
transcripts and avoid the introduction of artificial struc- 
tural changes, we used a method recently described by 
Holmquist et al. (43), based on the work of the 
Theobald-Dietrich and co-workers (44). Basically, we con- 
structed and synthesized a DNA template of sRNAi 62 for 
T7 polymerase, where the promotor was fused to a ham- 
merhead ribozyme (Supplementary Figure S3), further 
linked to the sequence of choice, here sRNA 162 . As soon 
as polymerization of the transcript starts, the RNA is sim- 
ultaneously folded into its hammerhead shape. When 
folding is completed, the ribozyme is autocatalytically 
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self-cleaved, whereas the transcription reaction of the 
downstream sequence continuous. Thus, the fusion guar- 
antees transcription at the native +1 site of sRNAi 62 , and 
can be engineered for any template choice, thereby 
providing an efficient method for the preparation of 
native RNA transcripts (Supplementary Figure S3). 

In vitro structure analysis 

Structure probing of sRNAs and mRNAs were conducted 
in total volume of 10 ul. Before RNase Tl and A cleavage 
the RNA (~0.1pmol) was denatured for lmin at 95°C 
and chilled on ice for 5min, then lOx Structure Buffer 
(Ambion; Applied Biosystems), 1 ug yeast RNA 
(Ambion; Applied Biosystems) were added. Following a 
renaturation step for lOmin at 37°C, 25 mM 
lead(II)acetate (Carl Roth, Karlsruhe), 0.01 U RNase Tl 
(Ambion; Applied Biosystems) or RNAse A (Ambion; 
Applied Biosystems) was added and incubated for to 2, 
3 and 5 min. The reactions were either stopped with 0.2 M 
EDTA and precipitated, or by directly adding gel loading 
buffer II (Ambion; Applied Biosystems). 

For RNase Tl ladders, the labeled RNA (~0.2pmol) 
was incubated for 1 min at 95°C in lx Sequencing Buffer 
(supplied with RNase Tl). Subsequently, 0.1 U RNase Tl 
was added to mix and further incubated for 5-10 min at 
37°C. OH ladders were generated by incubating ~0.2pmol 
labeled RNA for 5 min at 95°C in lx alkaline hydrolysis 
buffer (supplied with RNase Tl). Reactions were stopped 
as described earlier. The samples were analysed on 8% 
polyacryl amide/7M urea sequencing gels and visualized 
with a phosphoimager (FLA-5000 Series, Fuji). 

In vitro binding assays 

Electro mobility shift assays were conducted in a total 
volume of 10 ul in the presence of IX structure buffer 
(Ambion; Applied Biosystems) and 1 ug yeast RNA 
(Ambion; Applied Biosystems). 20pmol of in vitro 
transcribed RNA were dephosphorylated and radio- 
actively 5' labeled as described earlier. 5nM of the 
labeled RNA were incubated in presence with increasing 
amounts of the target RNA (8, 16, 32, 125, 250, 500, 1000, 
2000 nM) for 15 min at 37°C and subsequently separated 
on native 6-8% poly acrylamide gel in a 0.5x Tris-borate 
buffer system (0.45 M, pH 8.0). Gels were analysed using a 
phosphoimager (FLA-5000 Series, Fuji). 

Transcriptome analyses 

For genome-wide expression profiling genome wide micro- 
arrays representing 97% of the ORFs were used (41,45). 
Total RNA was extracted from M. mazei sRNA 162 - 
overexpressing cultures and the wt containing pWM321 
grown with 150mM methanol as carbon and energy 
source as described in Veit et al. (41). In general, RNA 
was extracted at turbidities OD 60 o = 0.5. Purified RNA 
was converted to cDNA and labeled by fluorescent Cy-3 
and Cy-5 dyes using the CyScribe first-strand cDNA 
labeling kit (GE Healthcare, Freiburg, Germany) as 
described (41,45). Microarray slides were incubated with 
aliquots of the cDNA preparations at 42° C overnight 
(Lucidea SlidePro hybridization chamber, GE 



Healthcare), for details of the hybridization and wash pro- 
cedures see Hovey et al. (45). Signal intensities were 
analysed using a GenePix 4 100 A scanner and data nor- 
malization and evaluation was performed using the 
GenePix Pro software version 6.0 (Axon Instruments, 
Union City, USA) as described recently (41,42). RNA 
isolated from three independent cultures were used inde- 
pendently in pairs of labeling reactions and, for one pair, a 
dye-swap experiment was also performed. Only if a tran- 
script had an abundance difference of at least 3-fold in the 
comparisons of wt RNA was the difference reported here 
to be significant. 

Quantitative reverse transcriptase PCR (qRT-PCR) 
assays were performed with a QuantiTect Probe 
RT-PCR Kit (Qiagen, Hilden, Germany) using a 7300 
real-time PCR system (ABI, Foster City, USA) and at 
least three independent RNA preparations for each 
strain as described (42,46). Primer sets used including 
the control genes are listed in the Supplementary Table 
S2. The fold change in transcript abundance for genes of 
interest was determined by comparison with the threshold 
cycle (C t ) of transcripts of three control genes (MM 1621, 
MM2181, MM1215). The fold change in the abundance of 
a transcript was calculated using the formula fold 
change = 2" AAO as described (47). 

Computational target predictions 

The ribosome binding site (RBS) sequence positions for all 
annotated M. mazei genes were predicted following the 
approach of (48). Significant RBS locations were 
obtained for 55.2% of the genes. Putative base pair inter- 
actions with sRNA 162 were predicted for all ORFs with at 
least 3-fold change in transcript level in the microarray 
analysis. For ORFs organized in an operon structure, 
the first ORF of the respective operon was also 
included. Sequences used for the interaction prediction 
contained the full length 5'-UTR (36) and additional 
100 nt of the CDS. If the exact 5'-UTR start was 
unknown, 200 nt upstream of the annotated translation 
start were used instead. Interaction predictions were 
computed with the tool IntaRNA (49,50) requiring an 
interaction seed of eight consecutive base pairs. All 
mRNA subsequences were folded locally in a 100 nt 
window with 50 nt maximal base pair span. For those 
ORFs which showed reduced transcript levels in the 
sRNAi 6 2 overproducing mutant using the microarray, 
we selected all interactions that were predicted to be 
located in the mRNA region from —39 in the 5'UTR to 
+ 19 in the CDS, which is the maximal region covered by 
ribosomes (51). For all ORFs with increased transcript 
levels, we studied the influence of the putative interaction 
on the accessibility of the predicted RBS sequence. All 
interactions predicted to be located upstream of the RBS 
sequence were included. A measure for the accessibility of 
a subsequence is its probability to be unpaired, i.e. free of 
intra-molecular base pairings. We computed the unpaired 
probability of the RBS sequence before and after the 
putative interaction of the mRNA with the sRNA 162 
(49). When the change in this unpaired probability was 
greater than 0.001, we assumed that the structural 
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rearrangement in the 5'-UTR caused by the interaction 
increases the accessibility of the RBS to activate the trans- 
lation of the gene. 

RESULTS 

Characterizing sRNAi 62 

We have recently identified sRNA 162 in our dRNA-seq 
approach to study the transcriptome of M. mazei (36). 
Because northern blot analysis revealed growth phase-de- 
pendent expression and processing of this particular 
sRNA, we selected sRNA 162 for further functional char- 
acterization and analysis of its target genes. The sRNA 162 
gene is located in the 739 bp intergenic region between 
MM2441, encoding a transcriptional regulator of the 
ArsR family, and MM2442, encoding a conserved 
protein of unknown function (Figure 1A). The TSS and 
termination site of sRNA 162 were determined by 5' RLM 
and 3' RACE, respectively (see 'Materials and Methods' 
section), and revealed a primary transcript of 191 nucleo- 
tides (nt) (Figure IB and C). In addition, the RACE 
approach revealed the presence of the homologous 
sRNA 171 with exactly the same transcript length, which 
is encoded in the intergenic region between MM2448 and 
MM2449 (discussed later). Expression of both RNAs was 
confirmed by northern blot analysis, and further revealed 
that the constitutively transcribed full length sRNAi 62 is 
processed into three different stable 5' fragments (nt 1-65, 
1-60, 1-55). In stationary growth phase, complete pro- 
cessing of sRNA 162 is observed independently of the 
supplied carbon source (Figure 1C; Supplementary 
Figure SI A). This is possibly achieved by exonucleases 
which degrade the 3' fragment, since the corresponding, 
processed 3' fragment(s) are neither detectable by northern 
blot analysis nor qRT-PCR (data not shown). In addition 
to the primary sRNA 162 transcript and its respective 5' 
processing fragments, other fragments were also detected 
using the probe directed against sRNA 162 as evident in a 
sRNA 162 deletion strain (Figure 1C). These are probably 
due to unspecific cross-hybridization (marked *) or 
cross-hybridization with sRNA 171 . To discriminate 
between the two homologous sRNAs, we additionally 
generated riboprobes specific for the respective full 
length sRNAs by in vitro transcription resulting in 
neglectable cross-hybridization in Northern blots 
(Supplementary Figure SIB and C). Hybridization with 
these riboprobes demonstrated that sRNAni is also con- 
stitutively expressed and appears to be immediately pro- 
cessed into stable 5' fragments, especially when cells enter 
stationary phase. However, additional processing frag- 
ments were detectable of ~150 and 120nt in length 
(Supplementary Figure SIC), indicating different process- 
ing of the two RNA species and/or different stabilities of 
the resulting 3' fragments, though high conservation of 
their primary sequence (Figure ID). 

BLAST searches for sRNA 162 homologs in other pro- 
karyotes detected the presence of homologous sequences 
exclusively in the two close relatives, M. acetivorans and 
M. barkeri. Secondary structure alignments of all 
sRNA 162 homologs revealed not only highly conserved 



structural elements [stem loop (SL) 1-3] but also a 39 nt 
single stranded linker region (SLR) between SL2 and SL3 
(Figure ID), representing the most variable part within the 
predicted structural alignment. In vitro secondary struc- 
ture probing using RNase Tl, RNase A- or Pbll- 
cleavage of in vitro synthesized 5'-end-labeled sRNA| 62 , 
confirmed the in silico predicted secondary structure with 
only marginal variations in the organization of internal 
loops of SL3 (compare Figures IE and 2). Analysis of 
the homologous sRNA 171 revealed a similar secondary 
structure consisting of three stem loops (SL1-SL3), with 
the two structured elements SL2 and SL3 interrupted by a 
SLR of 39 nt (Supplementary Figure S2). Single nucleo- 
tide exchang es m comparison with sRNA^6 2 are mostly 
present in the 5' part of sRNA 171 ranging from position 
1 to 80. Of those located in the structured parts, the 
majority represents compensatory base-pair exchanges or 
nucleotide exchanges within bulge loops which do not 
affect base-pairing [e.g. G42 (sRNA 162 ), A42 
(sRNA 17 i)]- Changes downstream of SL2 tend to be A 
to U substitutions, whereas at the cognate positions of 
sRNA 162 Gs and Cs are dominating (Figures ID and 2; 
Supplementary Figure S2). 

Identification of potential sRNAi 62 targets by genetic and 
computational approaches 

A deletion mutant of sRNAi 62 and a mutant expressing 
the sRNA 162 gene from the low copy plasmid pWM321 
under the control of its native promoter (further 
designated as 'over-expressing mutant') were generated 
(see 'Materials and Methods' section). Northern blot 
analysis confirmed deletion of sRNA 162 and showed 
~30-fold higher transcript levels in the over-expressing 
mutant strain (Figure 1C). To identify potential target 
genes of sRNA 162 , genome-wide changes in transcript 
levels in this mutant in comparison with the wt strain 
were studied using established genomic microarrays for 
M. mazei (41,45). The transcriptome analysis demons- 
trated that approximately 185 of the ORFs showed sig- 
nificantly different transcript levels (> 3-fold) in the 
over-expressing mutant compared with the wt strain. 
These include 48 genes involved in energy and construc- 
tional metabolism, two transcriptional regulators, 14 
genes encoding transport or membrane proteins and add- 
itional 51 genes with unknown function (Supplementary 
Table S3). Strikingly, elevated transcript levels of a high 
number of genes (7 operons) encoding soluble methyltran- 
sferases involved in degradation of methylamines were 
obtained in the mutant, although the strains were grown 
on methanol (Table 1). A transcriptome analysis of the 
sRNA 162 deletion mutant was not performed. Because 
the first 65 nt of sRNA 162 overlap with the 5'UTR of the 
flanking gene MM2442, we assumed polar effects in the 
deletion mutant on MM2442 (discussed later, second 
target). 

All ORFs with significant changes in transcript levels in 
the overexpression strain (>3-fold) were further analysed 
in silico to predict potential direct targets of sRNA 162 . The 
tool IntaRNA (49) was used to predict putative inter- 
actions between sRNA 162 and each of the respective 
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Figure 1. Characterization of sRNA 162 . (A) Genomic context of sRNA 162 . (B) Promotor and terminator region of sRNA 162 . The 5' end (+1) of 
sRNA 162 was determined by 5'-RLM RACE, as well as the 3' end (*) by 3'-RACE analysis. (C) Northern blot analysis of total RNA of wt and 
sRNA [62 mutant strains grown with methanol from early exponential phase (OD 600 of 0.2, lanes 1), mid exponential phase (OD 600 of 0.5, lanes 2) 
and stationary phase (OD 60 o of 0.7, lanes 3). The three 5'-fragments of sRNA I62 and the primary transcript are indicated by arrows. 
Cross-hybridization of the sRNA 162 probe is indicated by an asterisk. The lower panel shows the expression of 5S rRNA of the respective RNA 
preparation, wt, wild-type strain; sRNA 162 OE, sRNA 162 overexpressed from pWM321 in wt; and AsRNA| 62 , chromosomal sRNA 162 deletion 
mutant. (D) Secondary structure alignment of sRNA 162 homologues identified by computer-based searches of related Methanosarcina species per- 
formed with LocARNA (52), Ma, sRNA 162 homologue of M. acetivorans C2A; Mb, M. barkeri strain fusaro. (E) Consensus secondary structure 
prediction by RNalifold (53). SL, stem loop. 



mRNAs, including either their native 5'UTR, if known 
(36), or if not known, including 200 nt upstream of the pre- 
dicted translational start site (TLS) (see 'Materials and 
Methods' section). The obtained results for selected 
target candidates are summarized in Table 1. Additional 
potential interactions are listed in Supplementary 



Tables S4 and S5. The best scoring target was MM2441 
encoded in reverse orientation upstream of the sRNA 162 
gene (Figure 1A). The predicted interaction between the 
two RNAs is centered on the potential RBS and the TLS 
of mRNA2441 as depicted in Figures 3 A and B with a free 
energy score of — 17.6kcal/mol. This score is considerably 
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Table 1. Transcript levels of selected genes involved in TMA degradation and MM2441 of M. mazei sRNA 162 -overexpressing mutant versus 
M. mazei wt during growth on methanol as sole carbon and energy source identified by global expression profiling using genomic microarrays 
(41,45) 



Orf ID" Gene / Protein designation" 



Fold regulation 
(mutantversus wt) 



Interaction 
energy [kcal/mol] 



MM1435 Monomethylamine permease, mtmP\ 

MM1436 Monomethylamine:corrinoid methyltransferase, mtmB! (C-terminal domain) 

MM1437 Monomethylaminexorrinoid methyltransferase, mtmB\ (N-terminal domain) 

MM1438 Monomethylaminexorrinoid protein, mtmQ\ 

MM 1687 Dimethylamine:corrinoid protein, mtbCi 

MM1688 Trimethylaminexorrinoid methyltransferase, m»B| (N-terminal domain) 

MM 1689 Trimethylaminexorrinoid methyltransferase, mttB\ (C-terminal domain) 

MM1690 Trimethylaminexorrinoid protein, mttC\ 

MM1691 Trimethylamine permease, mtt~P\ 

MM1693 Dimethylaminexorrinoid methyltransferase, mtbB { (N-terminal domain) 

MM1694 Dimethylaminexorrinoid methyltransferase, mtbB\ 

MM2047 Trimethylaminexorrinoid protein, mttC 2 

MM2048 Trimethylaminexorrinoid methyltransferase, mttB 2 (N-terminal domain) 

MM2049 Trimethylaminexorrinoid methyltransferase, mttB 2 (C-terminal domain) 

MM2051 Dimethylaminexorrinoid methyltransferase, mtbB 2 (N-terminal domain) 

MM2052 Dimethylaminexorrinoid protein, mtbC 2 

MM2961 Dimethylaminexorrinoid protein, mtbCi 

MM2962 Dimethylaminexorrinoid methyltransferase, mtbB } (N-terminal domain) 

MM2963 Dimethylaminexorrinoid methyltransferase, mtbB 3 (C-terminal domain) 

MM3334 Monomethylaminexorrinoid protein, mtmC 2 

MM3335 Monomethylaminexorrinoid methyltransferase, mtmB 2 (N-terminal domain) 

MM3336 Monomethylaminexorrinoid methyltransferase, mtmB 2 (C-terminal domain) 

MM2441 Transcriptional regulator, ArsR family 



operon 



operon 



operon 



operon 



operon 



operon 



operon 



0.36 ± 0.01 
8.51 ± 0.94 

5.41 ± 0.71 
13.37 ± 1.91 
4.45 ± 0.57 
46.54 ± 5.19 
8.78 ± 3.57 

8.42 ± 3.09 
8.93 ± 4.56 

18.58 ± 0.81 
20.63 ± 4.76 

5.48 ± 2.85 
3.01 ± 1.48 
9.27 ± 1.97 

11.57 ± 3.66 
23.89 ± 15.25 

17.22 ± 1,06 
15.28 ± 2.13 
10.42 ± 4.88 

4.64 ± 2.29 
10.46 ± 3.50 
8.71 ± 2.38 

0.24 ± 0.01 



-6.58 
-8.04 



-1.11 



-2.67 



-3.73 



-17.56 



The in silico predicted interaction energy with sRNA|62 is given. 
! 'As defined in (54). 

b Gene induction is represented as the ratio of median. Mean values obtained from five microarray experiments that satisfy the criteria defined in the 
'Material and Methods' section are indicated. 



better than the second best ranked target (MM2339) with 
a free energy score of —10.9 kcal/mol (Supplementary 
Table S4), or than any of the other methyltransferase 
genes with altered expression upon sRNA 162 overe- 
xpression (see above and Table 1). Thus, the predicted 
primary target MM2441 was selected for further biochem- 
ical validation. 

In vitro verification of a direct sRNAi 62 /mRNA 2441 
interaction 

Transcript analysis by 5'-RLM RACE revealed that 
MM2441 is encoded as the second gene in a bicistronic 
operon together with MM2440 (encoding a hypothetical 
protein), which is preceded by a short 5'-UTR of 36 nt 
(Figure 3A). Furthermore, the two ORFs overlap 
because the predicted TLS of MM2441 is located within 
the 3' coding region of MM2440, 37 nt upstream of the 
translational stop of MM2440 (Figure 3A). Electro- 
phoretic mobility shift assays using a 5' labeled in vitro 
synthesized fragment of the bicistronic mRNA 
MM2440-41 (+1 to +469 relative to the TSS) including 
the predicted binding site (MM2440-41 LO ng; Figure 3A) 
and in vitro synthesized non-labeled full length sRNA 162 
(see 'Materials and Methods' section), demonstrated that 
binding to 5nmol MM2440-41 L ong occurred at concen- 
trations >32nM of sRNA 162 (Figure 3D). The reverse 



assay, using 5'-labeled sRNA 162 and non-labeled 
MM2440-41 L ong> further verified a direct interaction 
between the two RNAs (Supplementary Figure S4A). 
The interaction was strongly diminished when equal 
amounts of non-labeled sRNAi 62 were used as competitor 
RNA, confirming specificity of the binding (compare lines 
9 and 10, Supplementary Figure S4A). In contrast, up to 
2 uM of sRNA 162 did not affect the mobility of a shorter 5' 
fragment of the target mRNA, which did not include the 
predicted binding site (MM2440-41 

short? ranging 

from+1 to 248; Supplementary Figures S4B and C), 
demonstrating that the predicted interaction site is com- 
pulsory for binding. The absence of the predicted, 
non-structured binding site in sRNA 162 (sRNA 162A63 . 88 ) 
or changing the region into polyT (sRNA 162 M3), 
resulted in a complete loss of the mobility shift of 
MM2440-41 LO ng (Figures 3E and F; reverse assays 
Supplementary Figure S4D and E). Triple point mutations 
within the potential interacting site changing either 
position 64-66 in sRNA 162 from CAC antistart codon to 
GUG (Ml mutant) or position 79-81 from CAU to GUA 
(M2 mutant), further demonstrated that the M2 mutation 
totally abolished binding (Figure 3G; reverse assay Sup- 
plementary Figure S4F), whereas binding of the Ml 
mutant is only slightly reduced (Figure 3H and Supple- 
mentary Figure S4G reverse assay). Introducing 
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Figure 2. Structure mapping of 5'-end-labeled sRNA| 62 and proposed secondary structure of sRNAi 62 . (A) 5nM 5'-end-labeled sRNA^ was 
subjected to RNase Tl, lead(II) and RNase A cleavage. The cleavage was performed for 2, 3 or 5min, respectively. Lane OH: alkaline hydrolysis 
laddar. Lane Tl: RNase Tl laddar under denaturing conditions. The position of the cleaved Gs is given on the left of the gel. Lane C: untreated 
sRNA 162 . The approximate positions of stem-loop structures SL1, SL2 and SL3 according to the sRNA 162 structure shown in (B) are depicted by 
vertical bars on the right of the gel. (B) Proposed secondary structure of sRNA 162 determined by in vitro structure mapping. Cleavages according to 
(A) are indicated by circles for RNase Tl cleavage under both native and denaturing conditions, by squares for native RNase A cleavage and by 
black arrows for native lead(II) cleavage. The 3' ends of the three small processed fragments of sRNAi 62 detected by northern blot analysis are 
indicated as well as the predicted interacting site with MM2441 (black bar). Processing sites (PS1-3) are indicated. 



compensatory mutations at position 461 to 463 (M2') or 
446 to 448 (Ml') of MM2440-41 LO ng, totally restored the 
interaction with the M2 mutant sRNA, whereas the inter- 
action with sRNAi 6 2 Ml is significantly improved (Figure 
31 and J; Supplementary Figure S4H and I). These 
findings strongly indicate a specific binding of sRNA 162 
to the mRNA (MM2440-41), which most likely leads to 
masking of the RBS and the TLS of MM2441, as depicted 
in Figure 3B (wt). The binding affinity appears to be rela- 
tively low in comparison with a variety of bacterial trans- 
acting sRNAs (e.g. GcvB; (4)), indicating that a chaperone 
such as Hfq might be required to facilitate annealing of 
the two RNAs. Thus, the effect of Lsm-like proteins 
(MM339, MM2383) on binding has been tested by gel 
shift experiments in the presence of up to 75 uM purified 
his-tagged proteins (calculated for the hexamers); 
however, no significant effect was observed (data not 
shown). 



Gel shift assays using the 5' fragment (nt 1-65) 
of sRNA 162 (resembling the largest detectable stabile 
5' fragment in vivo) or the 3' fragment (nt 65-191) 
containing most of the binding site (nt 60-88) but 
missing the structured loops SL1 and 2 (Figure 3C) 
further demonstrated, that both of the truncated frag- 
ments are not able to interact with the target mRNA 
(Supplementary Figure S4J-M). Because both mutants 
lack parts of the binding site, this strongly suggests that 
the entire interacting site is essential for effective target 
binding, and/or the structured 5' part of sRNA 162 (SL1 
and 2), though not directly involved in the interaction, is 
highly required for a stable mRNA MM2440-41 target 
interaction in vivo. 

In vivo target validation by genetic approaches 

The mutant overexpressing sRNA 162 showed a partial 
growth defect when grown on methanol as sole carbon 
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Figure 3. sRNA| 6 2 interacts with the bicistronic mRNA MM2440-41 in vitro. (A) Overview of the transcriptional and translational features of 
MM2440-41. The boxed region indicates the predicted interacting site. Bold letters either represent the RBS or TLS of the respective ORFs. The 
in vitro transcribed target mRNAs are schematically drawn below. (B) In silico predictions of the sRNAi 62 -MM2440-41 interaction performed with 
IntaRNA (49). The wild-type (wt), triplet mutations (Ml and M2) and compensatory mutations (Ml' and M2') are indicated. (C) Cartoon of the 
proposed structure of sRNA 162 . Deletions or replacements of sRNA parts are indicated in red. (D-J) Electrophoretic mobility shift assays were 
performed using approximately 5nM of radioactively 5'-end-labeled MM2440-41 L ong (or compensatory mutants, respectively). The assays were 
performed as described in the 'Materials and Methods' section with increasing concentrations of unlabeled sRNA 162 or mutated species from 0 to 
2uM. After 15min incubation, samples were run on a native 6% PAA gel. The respective autoradiographs of the gels are shown. 



and energy source (Supplementary Figure S5). To gain a 
deeper insight, whether the growth reduction is due to 
interaction of sRNAi 62 with the bicistronic operon 
MM2440-41, or due to the interaction with a yet 
unknown target (or a combination of both), additional 
mutant strains were generated. Mutant derivatives of 
sRNA 162 identical to the variants used for EMSAs 
(Figure 3C), were ectopically expressed in the wt back- 
ground from pWM321 under the control of the native 
promoter and with the native terminator, including the 
separated 5' (1-65 nt) and 3' fragment (65-191 nt, contain- 
ing most of the predicted binding site), and the Ml to M3 
derivatives of sRNA 162 . Neither overexpression of the 
fragments nor the mutant derivatives of the full length 
sRNA 162 effected growth (Figure 4A); solely overexpres- 
sing full length wt sRNA 162 led to the observed significant 
change in growth rate. This finding strongly indicates that 
exclusively the primary transcript, containing the intact wt 
SLR, is capable to interact with the proposed target(s) and 
is crucial for the obtained growth defect, most likely due 
to down regulation of MM2441. However, at the current 
experimental status we cannot rule out the possibility that 
the phenotype is due to interactions with other potential 



targets, since the chromosomal compensatory mutations 
in MM2440-2441 nor a double deletion of sRNA 162 and 
MM2441 have not been generated due to experimental 
bottle necks (e.g. the lack of a second selection marker). 

To test whether processing of the sRNA is altered in 
case of the sRNA 162 derivatives northern blot analysis was 
performed (Figure 4B). Compared with the wt back- 
ground, overexpression and subsequent processing of 
full length sRNA 162 results in significant higher transcript 
levels of the primary transcript, as well as the three 5' 
fragments. A similar increase of the three 5' fragments 
was also detected upon overexpression of the 5' 
fragment (1-65 nt) of sRNA 162 (Figure 4B, lane 3) 
indicating that this extended 5' fragment is further 
processed into similar fragments and ratio as for the 
processing of the full length sRNA 162 . Remarkably, 
overexpression of the 3' fragment (nt 65-191), results in 
stable expression of the 3' portion without detectable pro- 
cessing pointing toward a crucial role of the structured 5' 
end of sRNA 162 for processing. However, the increased 
stability of the 3' fragment might also result from the arti- 
ficial truncation. Mutant derivatives of full length sRNA, 
Ml and M2, showed higher accumulation of the primary 
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time (h) 

Figure 4. Growth and northern blot analysis of sRNA 162 overexpressing mutants. (A) Growth of wt with additional sRNA 162 (sRNA 162 OE; empty 
diamond), 5' end of sRNA [62 (5' end; empty circles), 3' end of sRNAi 62 (3' end; filled squares), Ml mutant (Ml; empty squares), M2 mutant 
(M2; filled triangles) and M3 mutant (M3; filled circles), each ectopically overexpressed from pWM321 in wt background. Cells were grown on 
methanol to mid-exponential phase. (B) Northern blot analysis of the respective M. mazei strains using in vitro synthesized full-length RNA probes. 



transcript compared with the wt, strongly indicating 
altered processing. In contrast, when replacing the SLR 
by poly(T) (derivative M3), the primary transcript appears 
less stable compared with the Ml and M2 derivatives and 
the defined 5' fragments are still produced in high 
amounts, although as well some additional processing 
products are visible (Figure 4B). 

Considering that several soluble methyltransferases 
involved in degradation of methylamines showed 
elevated transcript levels upon sRNA 162 -overexpression 
(Table 1), growth was studied with trimethylamine 
(TMA) as sole energy and carbon (C) source. No differ- 
ence in expression or processing of sRNA 162 was obtained 
when cells were grown on TMA (Supplementary Figure 
SI A). Moreover, when grown on TMA the sRNA 162 
overexpression strain displayed the same growth pheno- 
type as observed with methanol (data not shown). 
However, when shifting cells from methanol to TMA, 
sRNA 162 -overexpression appears to result in significant 
faster adaptation following the C source shift compared 
with the wt (Figure 5) probably due to already synthesized 
soluble methyltransferases. To further support this sug- 
gestion, transcript levels of soluble methyltransferases in 
exponential growth phase with methanol were evaluated. 
Because of high sequence similarities of methyltransferase 
genes reliable transcript quantification using the PCR- 
fragment based microarray is not possible (45); thus, we 
performed qRT-PCR analysis using specific primers 
(46,55). A variety of methyltransferases with different sub- 
strate specificities showed significantly reduced mRNA 
levels of their respective operons, including methanol- 
dependent methyltransferases (nztaC.Bi, MMO 174-75; 
mtaC 3 B 3 , MM 1648-47), TMA-dependent methyltrans- 
ferases (m«Bid MM1688-90; mttB 2 C 2 , MM2049- 
47), dimethylamine-dependent methyltransferases 
(mtbdBu MM1687, 1693-94; mtbC 2 B 2 , MM2052-50) 
and monomethylamine-dependent methyltransferases 
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Figure 5. Growth behavior of the mutant overexpressing sRNA 162 
when shifted from methanol to trimethylamine in comparison with 
the wt. Cells were precultured as described in the 'Materials and 
Methods' section with methanol as sole carbon and energy source. 
After reaching exponential growth phase, cells were inoculated in 
fresh media either containing methanol (control) or trimethylamine 
(shift) as carbon and energy source. Growth was monitored by 
measuring the turbidity at 600 nm. Wt; sRNA| 62 OE, sRNA 162 
overexpressed from pWM321 in wt. Depicted are mean values and 
standard deviations of three biological replicates. 



(mfmCiB., MM1438-36) (Figure 6 and Supplementary 
Figure S6), explaining the reduced growth rates on 
methanol and TMA. Solely, mtmC 2 B 2 , encoding a 
monomethylamine-dependent methyltransferase 
(MM3334-36), showed approximately 5-fold increased 
transcript levels upon sRNAi 62 -overexpression, most 
likely leading to the observed faster adaptation of the 
sRNA 162 overproducing mutant strain after a shift from 
methanol to TMA (Figure 6). 

On the basis of these findings, we hypothesize that 
MM2441 is most likely post-transcriptionally regulated 
by sRNAi 62 and encodes a transcriptional regulator, 
which represses transcription of mtmB 2 C 2 encoding an es- 
sential soluble methyltransferase involved in TMA deg- 
radation. However, transcription of several other 
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Figure 6. Transcriptional profile of selected genes in exponential growth phase. For selected genes, the transcriptional levels were determined in 
exponential growth phase by qRT-PCR (for primers see Supplementary Table S5). Fold changes in the sRNA 162 overexpressing mutant (vs. wt) are 
given by mean values and standard deviation of three biologically independent experiments. MM0174, mtaCi encoding methanol corrinoid protein; 
MM1073, mtaC 2 encoding methanol corrinoid protein; MM1070, mtaA^ encoding methylcobalamin-coenzyme M methyltransferase; MM1438, 
mtmCi encoding monomethylamine corrinoid protein; MM1687, mtbC x encoding dimethylamine corrinoid protein; MM1690, mttCi encoding 
dimethylamine corrinoid protein; MM2047, mttCi encoding trimethylamine corrinoid protein; MM2052, mtbC 2 encoding dimethylamine corrinoid 
protein; MM2961, mtbC 3 encoding dimethylamine corrinoid protein; MM3334, monomethylamine corrinoid protein; MM2440, hypothetical protein; 
MM2441,/ transcriptional regulator, ArsR family; MM2442, hypothetical protein; MM2446, conserved protein. 



methyltransferase genes appears to be activated by 
MM2441 (see discussion). 

MM2442 as the second potential target of sRNAi 62 

Although short UTRs or even leaderless mRNAs have 
been predominantly found in (hyper-) thermophilic 
archaea (35,56), in methanoarchaea long 5' UTRs with 
an average size of ~150nt are more the rule than the 
exception (36). Because sRNAi 62 is located 150 nt 
upstream of the TLS of MM2442 in opposite orientation 
(Figure 1A), we determined the TSS of MM2442 encoding 
a conserved protein of unknown function. 5'-RLM RACE 
analysis identified the TSS of MM2442 to be 184nt 
upstream of the predicted TLS resulting in an overlap of 
65 nt in antisense orientation with the 5' end of sRNA 162 
(Figure 7A). This strongly suggests that besides targeting 
mRNA2440-41 in trans, sRNA 162 very likely represents 
also a c/.v-encoded antisense regulator of MM2442 
mRNA. Gel mobility shift experiments with full length 
sRNA 162 and its 3' and 5' fragments demonstrated that 
full length sRNA 162 indeed binds with its 5' region to the 
5'-UTR of MM2442 mRNA, whereas the isolated stable 5' 
fragment appears to have even a significant higher binding 
affinity (Figure 7B-D; reverse assays Supplementary 
Figure S7). Moreover, in vitro foot printing analysis con- 
firmed that sRNA 162 binds to the leader of MM2442 by its 
first 65 nt as depicted in Supplementary Figure S8. 



Characterizing sRNA m and its potential targets 

The homologous sRNA 17 i is encoded in the IGR between 
MM2448 encoding a transposase and MM2449, encoding 
a protein of unknown function, which however shows 
some homology to MM2442. Upstream of MM2448 a 
second transposase gene (MM2447) is present, which is 
further upstream flanked by MM2446, encoding a tran- 
scriptional regulator of the ArsR family with MM2441 as 
the closest homolog (Supplementary Figure S9). Overall, 
the high conservation of the two sRNAs (including their 
promoters) as well as the conservation of the flanking 
regions of sRNA 162 and sRNA m and their respective 
genomic organization suggests that sRNA 171 has been 
generated by a duplication of the complete gene locus 
MM2440-MM2442 most likely via a transposition event 
(Supplementary Figure S10). This is further supported by 
the finding that in the close relatives M. acetivorans and 
M. barkeri only a single copy of the sRNA is present 
in similarly organized genetic loci (Supplementary 
Figure S10). 

Even though sRNA 171 generally showed high sequence 
conservation to sRNAi 62 4 nt changes are obvious within 
the SLR, which might alter the binding properties of the 
sRNA (Figure ID). However, due to the similar genomic 
organization potential cross interaction between sRNA 171 
and the respective target mRNAs of sRNA 162 and vice 
versa were investigated by in silico predictions 
(Supplementary Figure S9C) and gel shift experiments. 
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Figure 7. The 5'-UTR of MM2442 as a sRNA 162 target. (A) Promotor region of MM2442 and sRNA 162 . The 5' end (+1) of MM2442 was 
determined by 5'-RLM RACE. The TATA box and BRE are indicated for both, MM2442 and the c/'s-encoded antisense sRNA| 62 . The predicted 
interacting site is boxed. (B-D) Electrophoretic mobility shift assays were performed using ~5nM of radioactively 5'-end-labeled MM2442. The 
reactions were performed as described in the 'Materials and Methods' section with increasing concentrations of unlabeled sRNA 162 or mutated 
species from 0 to 2uM. After 15min incubation, samples were run on a native 6% PAA gel. The respective autoradiographs of the gels are shown. 



This demonstrated a significantly lower binding affinity of 
sRNA 171 to the 2440-41 mRNA than observed for 
sRNA 162 (Supplementary Figure S9D; reverse assay 
Supplementary Figure S7) and no binding to MM2446 
for both sRNAs (Supplementary Figure S9F and G; 
reverse assays Supplementary Figure S7). Most interest- 
ingly, the second target of sRNA 162 , the 5' UTR of 
MM2442, was bound by sRNA 171 (Supplementary 
Figure S9E; reverse assay Supplementary Figure S7) 
with similar affinity as observed for the 5' fragment of 
sRNA 162 (Figure 7D). Consequently, the targets of both 
sRNAs generated by duplication might formerly have 
been the same, but following the transposition event 
both sRNA-genes and targets might have picked up 
several different mutations leading to diversification of 
targets and functions. 



DISCUSSION 

In the last two decades, a continuously growing number of 
sRNAs has been identified in prokaryotes. However, most 
of the functional characterization of sRNAs has been 
carried out in enterobacteria. Although expression of 
sRNAs has also been verified in several archaeal species 
(32-36,57-63), neither a potential target mRNA nor a mo- 
lecular mechanism of target regulation or physiological 
role has been elucidated in archaea so far. Considering 
that cellular processes in archaea in general show many 
mechanistic features more similar to their eukaryotic than 
bacterial counterparts, e.g. transcription and translation 
machineries (64-68), 3' targeting of mRNAs similar to 
eukaryotic miRNAs has been considered to be very 
likely for archaeal sRNAs. In this study, however, we 
present the first detailed characterization of an archaeal 
sRNA interacting with the 5' UTR of two mRNAs, a 
cis- and a fra«,y-encoded target. 

An archaeal sRNA masks the ribosome-binding site of its 
target 

We provide several lines of evidence, that sRNAi 62 acts in 
trans on the neighboring operon MM2440-41 by internally 



binding to the bicistronic mRNA. The interacting site of 
sRNA 162 was narrowed down to the SLR between SL2 
and 3 (Figure 2) by in vitro binding assays and in vivo 
studies, clearly demonstrating that the SLR is crucial for 
functionality (Figures 3 and 4A). However, the observed 
binding affinity of sRNA 162 to its target MM2441 is com- 
parably low (for comparison see refs. 4, 69, 70). No 
further enhancement was obtained in the presence of the 
two heterologously expressed Lsm-like proteins of 
M. mazei, although we currently cannot exclude 
non-physiological folding or incorrect assembly of Lsm 
oligomers due to heterologous expression. The finding 
that sRNAi 62 binds to the predicted RBS as well as the 
TLS of MM2441, strongly argues that the interaction 
inhibits translation initiation most likely resulting in 
dis-coordinated operon expression. Because we do not 
see any changes at the mRNA level, regulation probably 
most likely occurs at the translational level. The discovery 
that sRNA 162 targets the RBS within the 5' region of 
MM2441 mRNA elucidated an unexpected regulatory 
mechanism of an archaeal sRNA, which was up to now 
exclusively described for bacterial sRNAs. 

Homologs of RNase E and RNase III that are frequently 
involved in sRNA-mediated target regulation in bacteria 
(71-75) are not yet described in the archaeal domain. 
However, it is tempting to speculate that the increased 
turnover of sRNAi 62 into the stable 5' fragments 
observed in stationary phase (Figure 1C and D) might be 
facilitated by orthologous archaeal RNases. This 
growth-phase dependent processing of sRNA 162 leads to 
the accumulation of high amounts of stable 5' fragments, 
indicating that processing may result in another (regula- 
tory) outcome. However, ectopical overexpression of 
both, the longest 5' fragment (65 nt) and the 3' fragment 
of sRNA 162 (including most of the predicted binding 
site) does not affect growth as observed upon over- 
expression of the full length sRNAi 62 (Figure 4A), 
strongly suggesting an essential (stabilizing) effect of the 
5' fragment on sRNAi 62 , which appears crucial for func- 
tionality. Furthermore, it remains to be shown whether the 
resulting 5' fragments of sRNA 162 have additional regula- 
tory functions. 
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MM2441 most likely participates in the metabolic switch 
from methanol to (tri)methylamine 

In vivo evidence was obtained that sRNA 162 is involved in 
the adaptation process in response to different C-sources 
(e.g. shift from methanol to trimethylamine). Several 
soluble methyltransferases involved in methanol and 
methylamine-driven methanogenesis have been demon- 
strated to be differentially expressed upon overexpression 
of sRNAi62 but were not predicted as direct targets by 
bioinformatic predictions (see Supplementary Tables 
S3-S5). Based on our finding that sRNA 162 most likely 
interferes with MM2441 translation, we hypothesize that 
the indirect down-stream effects on transcription of 
several soluble methyltransferase genes in response to 
sRNA 162 overexpression results from the reduced synthe- 
sis of the ArsR-type transcriptional regulator MM2441. 
We further propose that during growth on methanol 
MM2441 represses transcription of the MM3334-3335 
operon, encoding a monomethylamine-dependent methyl- 
transferase (mtmB 2 ) and the cognate corrinoid protein 
(mtmC 2 ). However, other methyltransferase operons [e.g. 
MM1438-34 (m&wCjBjPi) and MM1687, 1693-94 
(mtbCiBi)] are most likely activated by MM2441. In line 
with this, significantly elevated transcript levels for the 
mtmB 2 Q 2 operon and reduced transcript levels for 
mtmCi and mtbC\ were detected by qRT-PCR in the 
sRNA 162 overexpressing mutant (Figure 6), probable 
due to decreased amounts of MM2441. This change in 
methyltransferase expression patterns upon sRNA 162 
overexpression most likely results in faster growth adap- 
tation after a switch from methanol to TMA (Figures 6 
and Supplementary Figure S6). 

Soluble methyltransferase genes involved in the degrad- 
ation of methanol and methylamines are in general among 
the most highly regulated genes in methanoarchaea 
(46,55,76-78). Several studies clearly demonstrated 
growth phase and carbon- and energy source dependent 



regulation of the respective enzymes in M. acetivorans and 
M. mazei and several overlapping regulatory circuits 
(55,76-79). A number of putative transcriptional regula- 
tors have been discovered, which impact on expression of 
the corresponding genes (76). For some of those, a dual 
functional role has been predicted, i.e. they act as activa- 
tors for transcription of several methyltransferase genes as 
well as repressors for others (76), which is in agreement 
with our findings. 

Overall, sRNA 162 appears to be involved in regulating 
the metabolic switch from methanol to methylamines 
during the transition of exponential to stationary growth 
phase, when cellular proteins and methylamines are 
degraded, most likely by post-transcriptional regulation 
of MM2441. We propose the following regulation 
(Figure 8): During exponential growth with methanol as 
carbon source, sRNA 162 is constitutively transcribed and 
subsequently binds to mRNA MM2440-2441. Thereby, 
sRNA I62 sequesters the translation initiation region 
(RBS and TLS) of MM2441, and controls its translation 
rates to ensure low MM2441 protein levels that are still 
sufficient to repress the mtmB 2 C 2 operon expression. 
Upon binding to its target, sRNAi 62 might be simultan- 
eously cleaved by RNases. However, in stationary growth 
phase the turnover of the sRNA 162 primary transcript is 
reinforced by a yet unknown factor that reflects entering 
the stationary growth phase (or artificially by 
overexpressing sRNA 162 ). This results in full translational 
repression of MM2441 and consequently de-repression of 
mtmB 2 C 2 transcription as well as reduced transcription of 
other methyltransferase operons (e.g. mtmC\BiPi and 
mtbCiBi) due to the absence of MM2441. 

sRNAi 62 interacts with cis- and fra/is-encoded targets via 
two distinct domains 

Despite acting as a trans-encoded antisense (as)RNA tar- 
geting MM2441, we further demonstrated that sRNA 162 



stationary phase 



derepression 




Figure 8. Proposed model of the sRNA 162 regulatory network. (A) sRNA 162 is constitutively expressed during exponential growth with methanol as 
sole carbon and energy source. When sRNA 162 is transcribed, it subsequently binds to mRNA MM2440-2441, sequestering the RBS and TLS of 
MM2441, thereby controlling translation rates to ensure constantly low MM2441 levels. The low MM2441 protein levels are sufficient to repress the 
mtmBC 2 operon expression. When the turnover of the primary transcript of sRNA 162 is reinforced (e.g. in stationary phase or by overexpressing 
sRNA 162 ) by a yet unknown factor, this results in enhanced translational repression of MM2441. At the same time, transcription of mtmB 2 C 2 is 
de-repressed by reduced MM2441 protein levels. The 5'- UTR of MM2442 encodes two potential small oligopeptides (sORF, black)). The regulatory 
outcome of the interaction with sRNA 162 is unknown. 
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also represents a czs-acting asRNA. While the SLR 
domain is essential for the proposed interaction with 
MM2441, we provided evidence that sRNA 162 targets 
the 5'-UTR of MM2442 by its very 5' end (nt 1-65). 
Because the transcript level of MM2442 is slightly nega- 
tively affected by sRNAi 62 overexpression (qRT-PCR, 
Figure 6), translational repression and subsequent degrad- 
ation of the transcript appears most likely, but transcrip- 
tional inference or promoter occlusion by sRNAi 62 is also 
possible. Upstream of the TLS of MM2442 (encoding a 
protein of so far unknown function) no RBS is present. A 
more detailed inspection of the long 5'-UTR of MM2442, 
further demonstrated that the 5' UTR sequence contains 
two small (s)ORFs, potentially encoding an oligopeptide 
of 23 or 30 amino acids (Supplementary Figure SI 1). This 
might indicate a coupled translation of the sORF and 
MM2442, and consequently, sRNA 162 could act by target- 
ing an upstream sORF of the bicistronic mRNA. Since 
sRNA 162 is masking the preceding RBS and CDS of the 
putative oligopeptide, it could lead to simultaneous trans- 
lational repression of the sORF and MM2442. A similar 
mechanism has been described for RyhB sRNA in E. coli 
(5), which represses the translation of the ferric uptake 
regulator (Fur) by binding to the RBS of a translationally 
coupled small ORF (designated as 'uORF') upstream of 
the Fur CDS. Moreover, translational activation of down- 
stream genes linked to a sORFs is mechanistically also 
possible, as exemplified by PhrS sRNA in Pseudomonas 
aeruginosa (80). The accumulation of the sRNA 162 5'-frag- 
ments during stationary growth probably maintains or 
even enhances translational repression of MM2442 and 
the potential small peptide. However, ectopic overproduc- 
tion of the 5' fragment of sRNA 162 did not result in an 
obvious phenotype. Therefore, the role of the sORF and 
MM2442 and the regulation by sRNA 162 is still elusive 
and remains to be investigated. 

Together these findings strongly suggest that sRNA 162 
acts on cis- and /ra«.s-encoded targets. Regulation of 
trans-encoded target mRNAs by ri.s-encoded asRNAs 
has been only proposed lacking any experimental 
evidence. Very recently Sayed et al. (81) described two 
RNAs organized as a type I toxin-antitoxin system con- 
sisting of SprAl, which encodes a small cytolytic peptide, 
and SprAl AS the ri.v-encoded riboregulator of SprAl ex- 
pression. Although SprAlAS overlaps with the 3' end of 
SprAl RNA, translation inhibition by SprAl is achieved 
by a functional domain outside of the complementary 
target sequence by forming a pseudoknot (81). However, 
sRNA 162 , differs from this system as it is not encoded in 
the same genomic locus as its rraw.s-encoded target 
mRNA, thus represents the first ris-encoded asRNA, add- 
itionally regulating a second target in trans. 

In conclusion, we have demonstrated for the first time 
that an archaeal sRNA, sRNAi 62 , has a regulatory 
function in archaea by interacting with the translation ini- 
tiation region (RBS and / or TLS) of its rra/w-encoded 
target(s), indicating that this archaeal sRNA acts similar 
as its bacterial counterparts. Moreover, we obtained 
evidence that sRNA 162 also targets a second, cw-encoded 
target. These obtained insights into the archaeal sRNA 162 



and its interaction with its targets blur the paradigm of a 
border between cis- and trans- encoded sRNAs. 
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